Intra-Landau-Level Cyclotron Resonance in Bilayer Graphene 
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Interaction driven integer quantum Hall effects are anticipated in graphene bilayers because of 
the near-degeneracy of the eight Landau levels which appear near the neutral system Fermi level. 
We predict that an intra-Landau-level cyclotron resonance signal will appear at some odd-integer 
filling factors, accompanied by collective modes which are nearly gapless and have approximate k 3 ^ 2 
dispersion. We speculate on the possibility of unusual localization physics associated with these 
modes. 
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Introduction — Because the Zeeman spin-splitting in 
most two-dimensional electron systems (2DES's) is much 
smaller than the Landau level separation, the magnetic 
band spectrum usually consists of narrowly-spaced dou- 
blets. When one of these doublets is half- filled and dis- 
order is weak, Coulomb interaction physics leads to fer- 
romagnetism i. e. to spontaneous spin polarization in the 
absence of a Zeeman field [3, 0, Q . In some circumstances 
Q other approximate Landau level degeneracies occur, 
often associated with layer degrees of freedom. These can 
also lead to broken symmetries which induce quasiparti- 
cle gaps and hence interaction driven integer quantum 
Hall effects. The case of bilayer 2DES's is particularly 
interesting because the which layer degree of freedom 
doubles Landau level degeneracies and leads to exciton 
condensation [B|, Q at odd filling factors and to canted 
anti-ferromagnetic states at even filling factors. In 
this Letter, we address the still richer case of graphene 
bilayer 2DES's in which chiral bands lead to an additional 
degeneracy doubling § at the Fermi energy of a neutral 
system. Bilayer graphene's Landau level octet is already 
apparent in present experiments [9( through the 8x (e 2 /h) 
Hall conductivity jump between well formed plateaus at 
Landau level filling factors v = —4 and v = +4. We 
anticipate that when external magnetic fields are strong 
enough or disorder is weak enough interactions will 
drive quantum Hall effects at the octet's seven interme- 
diate integer filling factors. We predict that these quan- 
tum Hall ferromagnets (QHFs) will exhibit unusual intra- 
Landau-level cyclotron modes at odd filling factors, and 
that the collective mode excitations at these filling fac- 
tors are nearly gapless even when there is no continuous 
symmetry breaking. Because the conductivity has Drudc 
weight centered near zero-energy, we speculate that local- 
ization physics and quantum-Hall related transport phe- 
nomena will also be anomalous. 

Graphene Bilayer Landau Levels — When trigonal warp- 
ing 11| and Zeeman coupling are neglected, the low 



and the influence of an external potential difference Ay 
between the layers is captured by 



energy properties of Bernal stacked unbalanced bilayer 
graphene are determined by electron-electron interac- 
tions and by a band Hamiltonian |8| Ti. = Ho + TL ex t 
where 

1 / n „rt2 \ 
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In Eqs. ([I])-®, 7? = p + (e/c)A is the 2D kinetic mo- 
mentum, 7r — tt x + iity, the 2x2 matrices act on the 
pseudospin degree of freedom associated with the two low 
energy sites [8[ (the top and bottom layer sites without a 
near-neighbor in the opposite layer), v is the single-layer 
Dirac velocity, 71 ~ 0.4eV is the inter-layer hopping am- 
plitude, and the effective mass m — ji/2v 2 w 0.054m e . 
H describes both K (£ = 1) and K' (f = -1) valleys 
provided that we choose the pseudospin representation 
(A, B) for K and (B, A) for K'. 



Defining the usual raising and lowering Landau level 
ladder operators a\a with = (Ib/V^H)^, where 
l B = (hc/eB) 1 / 2 = 25.6/ v /(fl[Tesla])nm is the mag- 
netic length, zero-energy eigenstates of Ho can be iden- 
tified using the property that a 2 (j> n — for 2D orbitals 
with Landau level index n = 0, 1. In bilayer graphene 
the n — and n = 1 orbital Landau levels are mem- 
bers of the same octet. This peculiarity is behind most 
of the physics explored in this paper. Neutral bilayer 
graphene's Landau-level octet is the direct product of 
three S — 1/2 doublets: real spin and which-layer 12 1 
pseudospins (as in a normal bilayer), and the Landau- 
level pseudospin n = 0, 1 degree of freedom which is re- 
sponsible for new physics. Zeeman coupling produces real 
spin-splitting Az while Ay gives rise to layer- split ting 
as in normal bilayers, but also to a small splitting of 
the Landau-level pseudospin which plays a central role 
in the physics: All — = Tholl where hio = 

2h 2 v 2 /l 2 B j! = 2.14B[Tesla] meV. 

Octet Hunds Rules — The octet HF Hamiltonian [l3} con- 
tains single-particle pseudospin splitting fields and direct 
and exchange interaction contributions: 
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FIG. 1: (Color online) Filling factor dependence of the inte- 
ger filling factor HF theory occupied state ( spectrum of the 
bilayer graphene octet at Ay = 0). Energies of occupied (red 
- solid lines) and unoccupied (blue - dashed lines) are in units 
of (7r/2) 1 / 2 e 2 /eis. The Zeeman field Az value in these units 
is 0.023 at a magnetic field of 20T. Octet space fractional 
pseudospin polarizations offset for clarity: spin(red boxes), 
valley(green circles) and LL pseudospin(blue triangles). 



occupied and empty states at the odd integer filling fac- 
tors of primary interest justify our weak-coupling theory. 
The octet filling, proceeding in integer increments start- 
ing from filling factor v = —4, follows a Hunds rule be- 
havior: first maximize spin-polarization, then maximize 
layer-polarization to the greatest extent possible, then 
maximize Landau-level polarization to the extent allowed 
by the first two rules. For balanced bilayers the layer 
symmetric states (S) are filled before the layer antisym- 
metric states (AS). The first four states to be filled are 
(S,n = 0, t),(S,n = 1, T),(AS,ri = 0, T) and (AS,n = 1, T) 
in this order. This sequence is then repeated for the next 
four states with down (J.) spin. The Hunds rules imply 
that the Landau-level pseudospin is polarized at all odd 
integer filling factors between v = —4 and v = 4. The 
physics of this new type of pseudospin polarization is the 
main focus of this paper. An important distinction be- 
tween layer and Landau-level polarization is that the for- 
mer is associated with spontaneous inter-layer phase co- 
herence whenever a Landau level occupies both layers si- 
multaneously, whereas the latter polarization is driven by 
the Landau-level dependence of the microscopic Hamil- 
tonian. 



where n = 0, 1 are LL indices, r, a = t(b) are valley in- 
dices, a, /3 =t (I) are spin-indices, and £ T ( a ) = 1(— 1) 
for t(b) layer and f (!) spins respectively. In Eq. © 
p T = J2na Prraa ^ s the total electron density in layer 
r. The density-matrix = (ci 2<J pC niTa ) must be 

determined self consistently by occupying the lowest en- 
ergy eigenvectors of H.hf- The Hartree- field Eh captures 
the electrostatic contribution to the bilayer capacitance, 
Eh = (e 2 /sIb)(<1/21b), and the exchange fields capture 
fermion quantum-statistics: 

/j2 
j^v $ (p)F n2n ,(p)F nni (-p). (4) 

In Eq.([4|) v± are the symmetric and antisymmetric com- 
binations of same (s) and different (d) layer electron- 
electron interactions (v s — 2Tre 2 /eq Vd — v s e~ qd ), 
and the form factors (F 00 (q) = er^ 1 ^ 2 ^, F 10 (q) = 
(iq x + q y )l B )e- ( - qlB)2/ VV2 = [F Q1 (-q)]* and F n (q) - 
(1 - {ql B f /2)e^ qlB)2/A ) reflect the character of the two 
different quantum cyclotron orbits. 

The solution of the Hartree-Fock equations for bal- 
anced bilayers (Ay = 0) is summarized in Fig.[[T] us- 
ing a Zeeman field strength corresponding to B = 20T. 
The large gaps (~ (n/8) 1 ^ 2 in e 2 /els units) separating 



Octet quantum Hall fcrromagnets have an interesting 
and intricate dependence on the external potential Ay. 
Because the two-layers are close together, a small value 
of Ay is sufficient to change the character of the layer 
polarization from the XY spontaneous-coherence form, 
to an Ising polarization form in which one layer is 
occupied before the other. We find that for Ay larger 
than a critical value Ay, the layer filling proceeds by 
filling the top layer first. (For v = -3, A v = 0.082(0.31) 
meV at B = 20(50) Tesla.) As we explain later, this 
filling sequence has qualitative consequences for the 
odd-integer filling factor LL pseudospin polarized states. 

Landau-Level Pseudospin Dipoles — We now focus on the 
LL pseudospin fluctuations of a state with odd-integer 
filling factor, freezing spin and layer degrees of freedom. 
The collective excitation spectrum of graphene bilayer 
octets as a function of v and Ay will be described in 
full detail elsewhere |l4{ . Fluctuating LL spinors are lin- 
ear combinations of n — orbitals (even with respect 
to their cyclotron orbit center) and n = 1 orbitals (odd 
with respect to orbit center) , and therefore carry an elec- 
tric dipole proportional to the in-plane component of 
their pseudospin. Because dipole-dipole interactions are 
long-range, they play a dominant role in the QHF long- 
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wavelength effective action Q. We find that 

S[m] = J dt[ J d 2 q A- d t fh- E[m}] , (5) 

where the first term is the Berry-phase contribution (3I Il5j 
and for small fluctuations away from m z = 1 (full n = 
polarization) 

E ^ = i;l^[^-A) 2 + ^(ml + ml)]. (6) 

where All = All/(g 2 /cIb)- The mass terms in Eq.([6]) 
are due to the single-particle splitting between n = and 
n = 1 levels and the interaction term is due to electric- 
dipole interactions. The absence of interaction contribu- 
tions to the mass terms is a surprise, since the interac- 
tion is Landau- level pseudospin dependent. We address 
this point below. Because of the in-plane electric dipolcs 
associated with LL pseudospinors, the long-wavelength 
pseudo-spinwave collective mode dispersion is not ana- 
lytic: hu -> (A 2 LL + A LL e 2 q/e) 1 / 2 , and for A LL -> 



is proportional to q 3 ^ 2 when exchange interactions are 
included in the energy functional. The in-plane dipoles 
are also responsible for the intra-Landau-level cyclotron 
resonance discussed below. 

To explain the absence of interaction contributions to 
the mass terms and address shorter- wavelength fluctua- 
tions it is necessary to derive the action microscopically. 
It is convenient to temporarily restrict fluctuations to 
one space direction by considering Landau-gauge states 
in which the LL pseudospins at different guiding centers 
X fluctuate independently: 



M*]> = II>oxcJx + zixc\ x )\0), (7) 
x 



where the spinor components z nX satisfy the normaliza- 
tion constraint |z x| 2 + l z ix| 2 = 1- The corresponding 
imaginary-time action is 



S[z, z) = Sb+£ 



Xi 



dry 



+ £A LL zi X z lx > 



XX' 



M 



where Sb is the Berry's phase term and £ = (ip[z]\(Tl + 
7~Lint)\'4>[ z \) 1S the energy functional. In Eq. (8} the 
direct(ff) and exchange(_F) energy contributions depend 
on the LL pseudospin labels, 

Ki;Z 2 ( x ) = ^-J ^v q F nini ( q )F n2 n 3 (~q)e-^ x , 
FZZ 2 ( X ) = ^EV 9H ,A« 3 (q)^ 2 n 4 (-q). (9) 

q 

This action can be identified as the Schwinger boson fl5(| 
coherent state path integral representation of a model 
with pseudospins at each guiding center. We can intro- 
duce a bosonic creation operator a\ LX corresponding to 
z n x and let £[z, z] — > Ti[a\a\. 

To analyze fluctuations around the HF mean field 
state, we use the linear spin wave approximation 

aax — » 1 - 7:a x a x a 1>x ^a x . (10) 



Taking the continuum limit l/L y ^2 x = J dX/(2irlB), 
the action describing harmonic fluctuations can be writ- 
ten in Fourier space as S = So + SS where 




(11) 



with 

Z q = \M e =^ _J dp{l _ £) Mq l BP )e^+ZA LL , 

X q = ^e^-jd P p 2 ( q l BP )e^, (12) 

In Eq.( [TTj) we have restored [l6| two-dimensional 
wavevectors to recognize the system's spatial anisotropy. 
The first and second terms in the above expressions cap- 
ture the direct(-ff) and exchange^) contributions re- 
spectively and Jo and J2 are the zeroth and second order 
Bessel functions. The quadratic action in Eq. (jTTJ) has 
the familiar Bogoliubov form and the energy dispersion 
of the collective mode is given by: 




As shown in Fig. this collective mode has a roton 
minimum at qIb ~ 2.3 and approaches the Hartree-Fock 
theory band splitting for q — > 00 as expected. [H The sur- 
; prising absence of interaction contributions to the gap at 
q = can be understood by examining the dependence of 
the uniform state interaction energy on global rotations 
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Ay = 20meV 
Ay = lOmeV 
A v = 0.082meV 
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FIG. 2: Collective mode u> q of the Landau-level pseudospin 
polarized state in units of interaction strength e 2 /ds = 
11.2-y/ B[Tesla] meV as a function of qls at different values 
of the external potential difference Ay at a magnetic field of 
20 T. The black(solid) line indicates the qIb — * oo asymptote 
for As = 0. 



in LL pseudospin space: 

28[z\_ _e^ 
el B 



2ko| 2 ki| 



(14) 



The factor in square brackets above is 1 — |zi| 4 /4, inde- 
pendent of z\ to quadratic order. Notice that because 
All < for v = — 1, 3 the absence of interaction contri- 
butions to the gap implies that the fully spin-polarized 
state is unstable. The ground state at these filling fac- 
tors is instead 14j an XY state with spontaneous phase 
order. 

Intra- Landau- Level Cyclotron Resonance — Finally we 
show that the octet QHF will exhibit unusual intra-LL 
cyclotron modes at odd filling factors, focusing on the 
fully polarized v — —3, 1 cases. The dynamical conduc- 
tivity a± — cr xx ± ia xy can be evaluated using linear 
response theory. The projection of the current operator, 
ji = dTL/d-Ki, onto the octet space can be expressed in 
terms of LL pseudospins: 



gA B , h 
m7i ^ v% 



Aext 

c l 



(*))> 



(15) 



where the ac electric field Ei — (l/c)dAf xt /dt. The ac 
conductivity (£ = 1) is most simply evaluated by solving 
the LL pseudospin equation of motion with the j ■ A ext 
coupling included in the energy functional. We find that 



1 



m7i i(ui ± loll) 



(16) 



In the absence of interactions the conductivity has intra- 
octet peaks at the LL band-splitting frequency loll, in 



addition to inter-Landau-level peaks which do not ap- 
pear in the projected theory. Trigonal warping is ex- 
pected |8[ to have little influence on electronic properties 
over the broad field range over which hl^ 1 > v^m. By 
performing explicit numerical calculations on the four- 
band model of bilayer graphene for typical bilayer quan- 
tum Hall parameters, a magnetic field strength of lOTesla 
and Ay ss lOmeV, we have verified that both the posi- 
tion and the oscillator strength of the intra-Landau level 
resonance are shifted by less than 2%. 

These low-frequency absorption peaks should be 
visible in microwave absorption experiments. The 
appearance of tunable low-frequency peaks in o~(lo) is a 
surprise that might be quite interesting from the point 
of view of the quantum Hall localization physics, even 
in systems for which disorder dominates interactions. 
In normal quantum Hall systems, peaks in a± appear 
near the characteristic inter-Landau-level energy oj c 
and the strong localization physics which leads to flat 
broad quantum Hall pleateaus occurs only in systems 
with uj c t > 1. We conjecture that one requirement for 
odd-integer filling factor plateaus within the graphene 
bilayer octet is that ujllt > 1- Since loll is proportional 
to Ay, the strength of the quantum Hall effect can be 
tuned by a gate voltage which doesn't influence either 
the system's disorder or its total carrier density. 
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